
clear all
cd "C:\Users\bgilbert_a\Dropbox\Supply Side Oil Gas\BigWellData"
set more off

use "DrillBitParityAnalysisData.dta", clear

* Table 3: Summary statistics for main variables
# delimit;
* Anadarko;
summ oil_targf1 oil_assocf1 oilf1 eprod3oilf1 etprod3oilf1 gas_targf1 gas_assoc_captf1 
		gasf1 eprod3gasf1 etprod3gasf1;
* Chautauqua Platform;
summ oil_targf5 oil_assocf5 oilf5 eprod3oilf5 etprod3oilf5 gas_targf5 gas_assoc_captf5 
		gasf5 eprod3gasf5 etprod3gasf5;
* East Texas;
summ oil_targf8 oil_assocf8 oilf8 eprod3oilf8 etprod3oilf8 gas_targf8 gas_assoc_captf8 
		gasf8 eprod3gasf8 etprod3gasf8;
* Fort Worth;
summ oil_targf9 oil_assocf9 oilf9 eprod3oilf9 etprod3oilf9 gas_targf9 gas_assoc_captf9 
		gasf9 eprod3gasf9 etprod3gasf9;
* Permian;
summ oil_targf14 oil_assocf14 oilf14 eprod3oilf14 etprod3oilf14 gas_targf14 gas_assoc_captf14 
		gasf14 eprod3gasf14 etprod3gasf14;
* Prices;
summ pwti_ce phh_ce;

# delimit;
* Table 4: IV Summary Statistics;
* Productivity variables are in logs. Need levels for summary statistics;
foreach var of varlist
		prod3gasfNan prod3gasfNch prod3gasfNet prod3gasfNfw prod3gasfNpr
		prod3oilfNan prod3oilfNch prod3oilfNet prod3oilfNfw prod3oilfNpr
		prod3gasf21  prod3oilf21  ln2_gas_pftb21  ln2_oil_pftb21
		totgaswell21 totoilwell21 totgas21 totoil21
		ln_gas_well21 ln_oil_well21 ln_gas_captf21 ln_oilf21 
		an2_gas ch2_gas et2_gas fw2_gas pr2_gas
		an2_oil ch2_oil et2_oil fw2_oil pr2_oil {;
		gen exp`var' = exp(`var');
		};

summ expprod3gasfNan-exptotoil21 brent cdd_dev hdd_dev sum12_cdd_dev sum12_hdd_dev;


* Regressions for Tables 6 and 7, as well as Appendix Tables;

*** Target initial; 

# delimit; 
set more off;
* Anadarko oil target initial;
ivreg2 ln_oil_targf1 L.ln_oil_targf1 L.prod3oilf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc  L.totoil21 L.totgaswell21 L.ln2_oil_pftb21 L.ln2_gas_pftb21  L.prod3oilfNan L.prod3gasf21 L.prod3gasfNan ), 
	first savefirst savefprefix(fs_ano) endog(lphh_ce lpwti_ce) orthog( L.totoil21 L.totgaswell21 L.ln2_oil_pftb21 L.ln2_gas_pftb21  L.prod3oilfNan L.prod3gasf21 L.prod3gasfNan ) bw(auto) robust;
est sto anop;
ivreg2 ln_oil_targf1 L.ln_oil_targf1 L.prod3oilf1 lphh_ce lpwti_ce, bw(auto) robust;
est sto anols; 
ivreg2 ln_oil_targf1 L.ln_oil_targf1 L.prod3oilf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto ano1;
ivreg2 ln_oil_targf1 L.ln_oil_targf1 L.prod3oilf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNan L.prod3gasf21 L.prod3gasfNan ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNan L.prod3gasf21 L.prod3gasfNan ) bw(auto) robust;
est sto ano2;
ivreg2 ln_oil_targf1 L.ln_oil_targf1 L.prod3oilf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto ano3;
ivreg2 ln_oil_targf1 L.ln_oil_targf1 L.prod3oilf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto ano4;
ivreg2 ln_oil_targf1 L.ln_oil_targf1 L.prod3oilf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.and_gas L.and_oil L.an2_gas L.an2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.and_gas L.and_oil L.an2_gas L.an2_oil) bw(auto) robust;
est sto ano5;

# delimit; 
set more off;
* Anadarko gas target initial;
ivreg2 ln_gas_targf1 L.ln_gas_targf1 L.prod3gasf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc  L.prod3gasfNan L.prod3oilfNan  L.totwell21), 
	first savefirst savefprefix(fs_ang) endog(lphh_ce lpwti_ce) orthog(L.prod3gasfNan L.prod3oilfNan  L.totwell21) bw(auto) robust;
est sto angp;
ivreg2 ln_gas_targf1 L.ln_gas_targf1 L.prod3gasf1 lphh_ce lpwti_ce, bw(auto) robust;
est sto angls; 
ivreg2 ln_gas_targf1 L.ln_gas_targf1 L.prod3gasf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto ang1;
ivreg2 ln_gas_targf1 L.ln_gas_targf1 L.prod3gasf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNan L.prod3gasf21 L.prod3gasfNan ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNan L.prod3gasf21 L.prod3gasfNan ) bw(auto) robust;
est sto ang2;
ivreg2 ln_gas_targf1 L.ln_gas_targf1 L.prod3gasf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto ang3;
ivreg2 ln_gas_targf1 L.ln_gas_targf1 L.prod3gasf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto ang4;
ivreg2 ln_gas_targf1 L.ln_gas_targf1 L.prod3gasf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.and_gas L.and_oil L.an2_gas L.an2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.and_gas L.and_oil L.an2_gas L.an2_oil) bw(auto) robust;
est sto ang5;

		

# delimit; 
set more off;
* Chautauqua oil target initial;
ivreg2 ln_oil_targf5 L.ln_oil_targf5 L.prod3oilf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc  L.totgas21 L.totoil21  L.prod3oilf21 L.prod3oilfNch L.prod3gasf21  ), 
	first savefirst savefprefix(fs_cho) endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21  L.prod3oilf21 L.prod3oilfNch L.prod3gasf21  ) bw(auto) robust;
est sto chop;
ivreg2 ln_oil_targf5 L.ln_oil_targf5 L.prod3oilf5 lphh_ce lpwti_ce, bw(auto) robust;
est sto chols; 
ivreg2 ln_oil_targf5 L.ln_oil_targf5 L.prod3oilf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto cho1;
ivreg2 ln_oil_targf5 L.ln_oil_targf5 L.prod3oilf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNch L.prod3gasf21 L.prod3gasfNch ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNch L.prod3gasf21 L.prod3gasfNch ) bw(auto) robust;
est sto cho2;
ivreg2 ln_oil_targf5 L.ln_oil_targf5 L.prod3oilf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto cho3;
ivreg2 ln_oil_targf5 L.ln_oil_targf5 L.prod3oilf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto cho4;
ivreg2 ln_oil_targf5 L.ln_oil_targf5 L.prod3oilf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.chq_gas L.chq_oil L.ch2_gas L.ch2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.chq_gas L.chq_oil L.ch2_gas L.ch2_oil) bw(auto) robust;
est sto cho5;
	
# delimit; 
set more off;
* Chautauqua gas target initial;
ivreg2 ln_gas_targf5 L.ln_gas_targf5 L.prod3gasf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.prod3oilf21 L.prod3gasfNch L.prod3oilfNch), 
	first savefirst savefprefix(fs_chg) endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.prod3oilf21 L.prod3gasfNch L.prod3oilfNch) bw(auto) robust;
est sto chgp;
ivreg2 ln_gas_targf5 L.ln_gas_targf5 L.prod3gasf5 lphh_ce lpwti_ce, bw(auto) robust;
est sto chgls; 
ivreg2 ln_gas_targf5 L.ln_gas_targf5 L.prod3gasf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto chg1;
ivreg2 ln_gas_targf5 L.ln_gas_targf5 L.prod3gasf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNch L.prod3gasf21 L.prod3gasfNch ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNch L.prod3gasf21 L.prod3gasfNch ) bw(auto) robust;
est sto chg2;
ivreg2 ln_gas_targf5 L.ln_gas_targf5 L.prod3gasf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto chg3;
ivreg2 ln_gas_targf5 L.ln_gas_targf5 L.prod3gasf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto chg4;
ivreg2 ln_gas_targf5 L.ln_gas_targf5 L.prod3gasf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.chq_gas L.chq_oil L.ch2_gas L.ch2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.chq_gas L.chq_oil L.ch2_gas L.ch2_oil) bw(auto) robust;
est sto chg5;



		
# delimit; 
set more off;
* East Texas oil target initial;
ivreg2 ln_oil_targf8 L.ln_oil_targf8 L.prod3oilf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc  L.prod3gasfNet L.totgas21), 
	first savefirst savefprefix(fs_eto) endog(lphh_ce lpwti_ce) orthog(L.prod3gasfNet L.totgas21 ) bw(auto) robust;
est sto etop;
ivreg2 ln_oil_targf8 L.ln_oil_targf8 L.prod3oilf8 lphh_ce lpwti_ce, bw(auto) robust;
est sto etols; 
ivreg2 ln_oil_targf8 L.ln_oil_targf8 L.prod3oilf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto eto1;
ivreg2 ln_oil_targf8 L.ln_oil_targf8 L.prod3oilf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNet L.prod3gasf21 L.prod3gasfNet ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNet L.prod3gasf21 L.prod3gasfNet ) bw(auto) robust;
est sto eto2;
ivreg2 ln_oil_targf8 L.ln_oil_targf8 L.prod3oilf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto eto3;
ivreg2 ln_oil_targf8 L.ln_oil_targf8 L.prod3oilf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto eto4;
ivreg2 ln_oil_targf8 L.ln_oil_targf8 L.prod3oilf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.etx_gas L.etx_oil L.et2_gas L.et2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.etx_gas L.etx_oil L.et2_gas L.et2_oil) bw(auto) robust;
est sto eto5;	
	
# delimit; 
set more off;
* East Texas gas target initial;
ivreg2 ln_gas_targf8 L.ln_gas_targf8 L.prod3gasf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3gasf21 L.totgas21 ), 
	first savefirst savefprefix(fs_etg) endog(lphh_ce lpwti_ce) orthog(L.prod3gasf21 L.totgas21 ) bw(auto) robust;
est sto etgp;
ivreg2 ln_gas_targf8 L.ln_gas_targf8 L.prod3gasf8 lphh_ce lpwti_ce, bw(auto) robust;
est sto etgls; 
ivreg2 ln_gas_targf8 L.ln_gas_targf8 L.prod3gasf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto etg1;
ivreg2 ln_gas_targf8 L.ln_gas_targf8 L.prod3gasf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNet L.prod3gasf21 L.prod3gasfNet ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNet L.prod3gasf21 L.prod3gasfNet ) bw(auto) robust;
est sto etg2;
ivreg2 ln_gas_targf8 L.ln_gas_targf8 L.prod3gasf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto etg3;
ivreg2 ln_gas_targf8 L.ln_gas_targf8 L.prod3gasf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto etg4;
ivreg2 ln_gas_targf8 L.ln_gas_targf8 L.prod3gasf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.etx_gas L.etx_oil L.et2_gas L.et2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.etx_gas L.etx_oil L.et2_gas L.et2_oil) bw(auto) robust;
est sto etg5;


					
# delimit; 
set more off;
* Fort Worth oil target initial;
ivreg2 ln_oil_targf9 L.ln_oil_targf9 L.prod3oilf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_oil_pftb21 L.prod3oilf21 L.prod3gasfNfw  ), 
	first savefirst savefprefix(fs_fwo) endog(lphh_ce lpwti_ce) orthog(L.ln2_oil_pftb21 L.prod3oilf21 L.prod3gasfNfw  ) bw(auto) robust;
est sto fwop;
ivreg2 ln_oil_targf9 L.ln_oil_targf9 L.prod3oilf9 lphh_ce lpwti_ce, bw(auto) robust;
est sto fwols; 
ivreg2 ln_oil_targf9 L.ln_oil_targf9 L.prod3oilf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto fwo1;
ivreg2 ln_oil_targf9 L.ln_oil_targf9 L.prod3oilf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNfw L.prod3gasf21 L.prod3gasfNfw ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNfw L.prod3gasf21 L.prod3gasfNfw ) bw(auto) robust;
est sto fwo2;
ivreg2 ln_oil_targf9 L.ln_oil_targf9 L.prod3oilf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto fwo3;
ivreg2 ln_oil_targf9 L.ln_oil_targf9 L.prod3oilf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto fwo4;
ivreg2 ln_oil_targf9 L.ln_oil_targf9 L.prod3oilf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.fwth_gas L.fwth_oil L.fw2_gas L.fw2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.fwth_gas L.fwth_oil L.fw2_gas L.fw2_oil) bw(auto) robust;
est sto fwo5;

# delimit; 
set more off;
* Fort Worth gas target initial;
ivreg2 ln_gas_targf9 L.ln_gas_targf9 L.prod3gasf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc   L.prod3gasfNfw  L.prod3gasf21 ), 
	first savefirst savefprefix(fs_fwg) endog(lphh_ce lpwti_ce) orthog(L.prod3gasfNfw  L.prod3gasf21 ) bw(auto) robust;
est sto fwgp;
ivreg2 ln_gas_targf9 L.ln_gas_targf9 L.prod3gasf9 lphh_ce lpwti_ce, bw(auto) robust;
est sto fwgls; 
ivreg2 ln_gas_targf9 L.ln_gas_targf9 L.prod3gasf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto fwg1;
ivreg2 ln_gas_targf9 L.ln_gas_targf9 L.prod3gasf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNfw L.prod3gasf21 L.prod3gasfNfw ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNfw L.prod3gasf21 L.prod3gasfNfw ) bw(auto) robust;
est sto fwg2;
ivreg2 ln_gas_targf9 L.ln_gas_targf9 L.prod3gasf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto fwg3;
ivreg2 ln_gas_targf9 L.ln_gas_targf9 L.prod3gasf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto fwg4;
ivreg2 ln_gas_targf9 L.ln_gas_targf9 L.prod3gasf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.fwth_gas L.fwth_oil L.fw2_gas L.fw2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.fwth_gas L.fwth_oil L.fw2_gas L.fw2_oil) bw(auto) robust;
est sto fwg5;
					
					
					
# delimit; 
set more off;
* Permian oil target initial;
ivreg2 ln_oil_targf14 L.ln_oil_targf14 L.prod3oilf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc  L.prod3gasfNpr L.prod3gasf21 ), 
	first savefirst savefprefix(fs_pro) endog(lphh_ce lpwti_ce) orthog(L.prod3gasf21 L.prod3gasfNpr  ) bw(auto) robust;
est sto prop;
ivreg2 ln_oil_targf14 L.ln_oil_targf14 L.prod3oilf14  lphh_ce lpwti_ce, bw(auto) robust;
est sto prols; 
ivreg2 ln_oil_targf14 L.ln_oil_targf14 L.prod3oilf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto pro1;
ivreg2 ln_oil_targf14 L.ln_oil_targf14 L.prod3oilf14  
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNpr L.prod3gasf21 L.prod3gasfNpr ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNpr L.prod3gasf21 L.prod3gasfNpr ) bw(auto) robust;
est sto pro2;
ivreg2 ln_oil_targf14 L.ln_oil_targf14 L.prod3oilf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto pro3;
ivreg2 ln_oil_targf14 L.ln_oil_targf14 L.prod3oilf14  
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto pro4;
ivreg2 ln_oil_targf14 L.ln_oil_targf14 L.prod3oilf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prm_gas L.prm_oil L.pr2_gas L.pr2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.prm_gas L.prm_oil L.pr2_gas L.pr2_oil) bw(auto) robust;
est sto pro5;
	
# delimit; 
set more off;
* Permian gas target initial;
ivreg2 ln_gas_targf14 L.ln_gas_targf14 L.prod3gasf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3gasfNpr), 
	first savefirst savefprefix(fs_prg) endog(lphh_ce lpwti_ce) orthog(L.prod3gasfNpr ) bw(auto) robust;
est sto prgp;
ivreg2 ln_gas_targf14 L.ln_gas_targf14 L.prod3gasf14 lphh_ce lpwti_ce, bw(auto) robust;
est sto prgls; 
ivreg2 ln_gas_targf14 L.ln_gas_targf14 L.prod3gasf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto prg1;
ivreg2 ln_gas_targf14 L.ln_gas_targf14 L.prod3gasf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNpr L.prod3gasf21 L.prod3gasfNpr ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNpr L.prod3gasf21 L.prod3gasfNpr ) bw(auto) robust;
est sto prg2;
ivreg2 ln_gas_targf14 L.ln_gas_targf14 L.prod3gasf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto prg3;
ivreg2 ln_gas_targf14 L.ln_gas_targf14 L.prod3gasf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto prg4;
ivreg2 ln_gas_targf14 L.ln_gas_targf14 L.prod3gasf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prm_gas L.prm_oil L.pr2_gas L.pr2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.prm_gas L.prm_oil L.pr2_gas L.pr2_oil) bw(auto) robust;
est sto prg5;	



* Now estimate sum of our 5 basins and all basins in DI;

* All five of our basins aggregated;
# delimit;
ivreg2 oil_targfot L.oil_targfot L.prod3oilfot
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd L.totgas21 L.prod3gasf21),
	first savefirst savefprefix(fs_5o) endog(lphh_ce lpwti_ce) orthog(lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd   ) bw(auto) robust;
est sto otoilp;
ivreg2 gas_targfot L.gas_targfot L.prod3gasfot
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd L.totoil21 L.totgas21  L.prod3gasf21),
	first savefirst savefprefix(fs_5g) endog(lphh_ce lpwti_ce) bw(auto) robust;
est sto otgasp;

* sum of all basins in DI; 
# delimit;
ivreg2 oil_targftot L.oil_targftot L.prod3oilftot
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd   ),
	first savefirst savefprefix(fs_allo) endog(lphh_ce lpwti_ce) orthog(lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd ) bw(auto) robust;
est sto alloilp;
ivreg2 gas_targftot L.gas_targftot L.prod3gasftot
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd  ),
	first savefirst savefprefix(fs_allg) endog(lphh_ce lpwti_ce) bw(auto) robust;
est sto allgasp;
	
					
/*
Now for total initial production;
*/

# delimit; 
set more off;
* Anadarko oil total initial;
ivreg2 ln_oilf1 L.ln_oilf1 L.tprod3oilf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totoil21 L.prod3gasfNan L.totgaswell21 ), 
	first savefirst savefprefix(fs_anot) endog(lphh_ce lpwti_ce) orthog( L.totoil21 L.prod3gasfNan L.totgaswell21 ) bw(auto) robust;
est sto anopt;
ivreg2 ln_oilf1 L.ln_oilf1 L.tprod3oilf1 lphh_ce lpwti_ce, bw(auto) robust;
est sto anolst; 
ivreg2 ln_oilf1 L.ln_oilf1 L.tprod3oilf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto ano1t;
ivreg2 ln_oilf1 L.ln_oilf1 L.tprod3oilf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNan L.prod3gasf21 L.prod3gasfNan ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNan L.prod3gasf21 L.prod3gasfNan ) bw(auto) robust;
est sto ano2t;
ivreg2 ln_oilf1 L.ln_oilf1 L.tprod3oilf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto ano3t;
ivreg2 ln_oilf1 L.ln_oilf1 L.tprod3oilf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto ano4t;
ivreg2 ln_oilf1 L.ln_oilf1 L.tprod3oilf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.and_gas L.and_oil L.an2_gas L.an2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.and_gas L.and_oil L.an2_gas L.an2_oil) bw(auto) robust;
est sto ano5t;

# delimit; 
set more off;
* Anadarko gas total initial;
ivreg2 ln_gas_captf1 L.ln_gas_captf1 L.tprod3gasf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent L.prod3gasfNan L.prod3oilfNan), 
	first savefirst savefprefix(fs_angt) endog(lphh_ce lpwti_ce) orthog(L.prod3gasfNan L.prod3oilfNan) bw(auto) robust;
est sto angpt;
ivreg2 ln_gas_captf1 L.ln_gas_captf1 L.tprod3gasf1 lphh_ce lpwti_ce, bw(auto) robust;
est sto anglst; 
ivreg2 ln_gas_captf1 L.ln_gas_captf1 L.tprod3gasf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto ang1t;
ivreg2 ln_gas_captf1 L.ln_gas_captf1 L.tprod3gasf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNan L.prod3gasf21 L.prod3gasfNan ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNan L.prod3gasf21 L.prod3gasfNan ) bw(auto) robust;
est sto ang2t;
ivreg2 ln_gas_captf1 L.ln_gas_captf1 L.tprod3gasf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto ang3t;
ivreg2 ln_gas_captf1 L.ln_gas_captf1 L.tprod3gasf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto ang4t;
ivreg2 ln_gas_captf1 L.ln_gas_captf1 L.tprod3gasf1 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.and_gas L.and_oil L.an2_gas L.an2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.and_gas L.and_oil L.an2_gas L.an2_oil) bw(auto) robust;
est sto ang5t;

		

# delimit; 
set more off;
* Chautauqua oil total initial;
ivreg2 ln_oilf5 L.ln_oilf5 L.tprod3oilf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgaswell21 L.prod3gasfNch), 
	first savefirst savefprefix(fs_chot) endog(lphh_ce lpwti_ce) orthog(L.totgaswell21 L.prod3gasfNch) bw(auto) robust;
est sto chopt;
ivreg2 ln_oilf5 L.ln_oilf5 L.tprod3oilf5 lphh_ce lpwti_ce, bw(auto) robust;
est sto cholst; 
ivreg2 ln_oilf5 L.ln_oilf5 L.tprod3oilf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto cho1t;
ivreg2 ln_oilf5 L.ln_oilf5 L.tprod3oilf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNch L.prod3gasf21 L.prod3gasfNch ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNch L.prod3gasf21 L.prod3gasfNch ) bw(auto) robust;
est sto cho2t;
ivreg2 ln_oilf5 L.ln_oilf5 L.tprod3oilf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto cho3t;
ivreg2 ln_oilf5 L.ln_oilf5 L.tprod3oilf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto cho4t;
ivreg2 ln_oilf5 L.ln_oilf5 L.tprod3oilf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.chq_gas L.chq_oil L.ch2_gas L.ch2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.chq_gas L.chq_oil L.ch2_gas L.ch2_oil) bw(auto) robust;
est sto cho5t;	
	
# delimit; 
set more off;
* Chautauqua gas total initial;
ivreg2 ln_gas_captf5 L.ln_gas_captf5 L.tprod3gasf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent  blhdd blcdd bsh bsc L.prod3gasfNch L.prod3oilfNch L.totgas21 ), 
	first savefirst savefprefix(fs_chgt) endog(lphh_ce lpwti_ce) orthog(L.prod3gasfNch L.prod3oilfNch L.totgas21) bw(auto) robust;
est sto chgpt;
ivreg2 ln_gas_captf5 L.ln_gas_captf5 L.tprod3gasf5 lphh_ce lpwti_ce, bw(auto) robust;
est sto chglst; 
ivreg2 ln_gas_captf5 L.ln_gas_captf5 L.tprod3gasf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto chg1t;
ivreg2 ln_gas_captf5 L.ln_gas_captf5 L.tprod3gasf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNch L.prod3gasf21 L.prod3gasfNch ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNch L.prod3gasf21 L.prod3gasfNch ) bw(auto) robust;
est sto chg2t;
ivreg2 ln_gas_captf5 L.ln_gas_captf5 L.tprod3gasf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto chg3t;
ivreg2 ln_gas_captf5 L.ln_gas_captf5 L.tprod3gasf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto chg4t;
ivreg2 ln_gas_captf5 L.ln_gas_captf5 L.tprod3gasf5 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.chq_gas L.chq_oil L.ch2_gas L.ch2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.chq_gas L.chq_oil L.ch2_gas L.ch2_oil) bw(auto) robust;
est sto chg5t;
		

		
# delimit; 
set more off;
* East Texas oil total initial;
ivreg2 ln_oilf8 L.ln_oilf8 L.tprod3oilf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3gasfNet L.ln2_gas_pftb21 L.totgaswell21  ), 
	first savefirst savefprefix(fs_etot) endog(lphh_ce lpwti_ce) orthog(L.prod3gasfNet L.ln2_gas_pftb21 L.totgaswell21  ) bw(auto) robust;
est sto etopt;
ivreg2 ln_oilf8 L.ln_oilf8 L.tprod3oilf8 lphh_ce lpwti_ce, bw(auto) robust;
est sto etolst; 
ivreg2 ln_oilf8 L.ln_oilf8 L.tprod3oilf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto eto1t;
ivreg2 ln_oilf8 L.ln_oilf8 L.tprod3oilf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNet L.prod3gasf21 L.prod3gasfNet ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNet L.prod3gasf21 L.prod3gasfNet ) bw(auto) robust;
est sto eto2t;
ivreg2 ln_oilf8 L.ln_oilf8 L.tprod3oilf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto eto3t;
ivreg2 ln_oilf8 L.ln_oilf8 L.tprod3oilf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto eto4t;
ivreg2 ln_oilf8 L.ln_oilf8 L.tprod3oilf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.etx_gas L.etx_oil L.et2_gas L.et2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.etx_gas L.etx_oil L.et2_gas L.et2_oil) bw(auto) robust;
est sto eto5t;	
	
# delimit; 
set more off;
* East Texas gas total initial;
ivreg2 ln_gas_captf8 L.ln_gas_captf8 L.tprod3gasf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3gasfNet L.prod3oilfNet L.totoil21), 
	first savefirst savefprefix(fs_etgt) endog(lphh_ce lpwti_ce) orthog(L.prod3gasfNet L.prod3oilfNet L.totoil21) bw(auto) robust;
est sto etgpt;
ivreg2 ln_gas_captf8 L.ln_gas_captf8 L.tprod3gasf8 lphh_ce lpwti_ce, bw(auto) robust;
est sto etglst; 
ivreg2 ln_gas_captf8 L.ln_gas_captf8 L.tprod3gasf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto etg1t;
ivreg2 ln_gas_captf8 L.ln_gas_captf8 L.tprod3gasf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNet L.prod3gasf21 L.prod3gasfNet ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNet L.prod3gasf21 L.prod3gasfNet ) bw(auto) robust;
est sto etg2t;
ivreg2 ln_gas_captf8 L.ln_gas_captf8 L.tprod3gasf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto etg3t;
ivreg2 ln_gas_captf8 L.ln_gas_captf8 L.tprod3gasf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto etg4t;
ivreg2 ln_gas_captf8 L.ln_gas_captf8 L.tprod3gasf8 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.etx_gas L.etx_oil L.et2_gas L.et2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.etx_gas L.etx_oil L.et2_gas L.et2_oil) bw(auto) robust;
est sto etg5t;


					
# delimit; 
set more off;
* Fort Worth oil total initial;
ivreg2 ln_oilf9 L.ln_oilf9 L.tprod3oilf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent L.prod3gasfNfw), 
	first savefirst savefprefix(fs_fwot) endog(lphh_ce lpwti_ce) orthog(L.prod3gasfNfw  ) bw(auto) robust;
est sto fwopt;
ivreg2 ln_oilf9 L.ln_oilf9 L.tprod3oilf9 lphh_ce lpwti_ce, bw(auto) robust;
est sto fwolst; 
ivreg2 ln_oilf9 L.ln_oilf9 L.tprod3oilf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto fwo1t;
ivreg2 ln_oilf9 L.ln_oilf9 L.tprod3oilf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNfw L.prod3gasf21 L.prod3gasfNfw ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNfw L.prod3gasf21 L.prod3gasfNfw ) bw(auto) robust;
est sto fwo2t;
ivreg2 ln_oilf9 L.ln_oilf9 L.tprod3oilf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto fwo3t;
ivreg2 ln_oilf9 L.ln_oilf9 L.tprod3oilf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto fwo4t;
ivreg2 ln_oilf9 L.ln_oilf9 L.tprod3oilf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.fwth_gas L.fwth_oil L.fw2_gas L.fw2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.fwth_gas L.fwth_oil L.fw2_gas L.fw2_oil) bw(auto) robust;
est sto fwo5t;
	
# delimit; 
set more off;
* Fort Worth gas total initial;
ivreg2 ln_gas_captf9 L.ln_gas_captf9 L.tprod3gasf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3gasfNfw L.prod3gasf21 L.totoil21 L.totgas21), 
	first savefirst savefprefix(fs_fwgt) endog(lphh_ce lpwti_ce) orthog(L.prod3gasfNfw L.prod3gasf21 L.totoil21 L.totgas21) bw(auto) robust;
est sto fwgpt;
ivreg2 ln_gas_captf9 L.ln_gas_captf9 L.tprod3gasf9 lphh_ce lpwti_ce, bw(auto) robust;
est sto fwglst; 
ivreg2 ln_gas_captf9 L.ln_gas_captf9 L.tprod3gasf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto fwg1t;
ivreg2 ln_gas_captf9 L.ln_gas_captf9 L.tprod3gasf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNfw L.prod3gasf21 L.prod3gasfNfw ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNfw L.prod3gasf21 L.prod3gasfNfw ) bw(auto) robust;
est sto fwg2t;
ivreg2 ln_gas_captf9 L.ln_gas_captf9 L.tprod3gasf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto fwg3t;
ivreg2 ln_gas_captf9 L.ln_gas_captf9 L.tprod3gasf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto fwg4t;
ivreg2 ln_gas_captf9 L.ln_gas_captf9 L.tprod3gasf9 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.fwth_gas L.fwth_oil L.fw2_gas L.fw2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.fwth_gas L.fwth_oil L.fw2_gas L.fw2_oil) bw(auto) robust;
est sto fwg5t;
					
					
					
# delimit; 
set more off;
* Permian oil total initial;
ivreg2 ln_oilf14 L.ln_oilf14 L.tprod3oilf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	first savefirst savefprefix(fs_prot) endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto propt;
ivreg2 ln_oilf14 L.ln_oilf14 L.tprod3oilf14  lphh_ce lpwti_ce, bw(auto) robust;
est sto prolst; 
ivreg2 ln_oilf14 L.ln_oilf14 L.tprod3oilf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto pro1t;
ivreg2 ln_oilf14 L.ln_oilf14 L.tprod3oilf14  
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNpr L.prod3gasf21 L.prod3gasfNpr ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNpr L.prod3gasf21 L.prod3gasfNpr ) bw(auto) robust;
est sto pro2t;
ivreg2 ln_oilf14 L.ln_oilf14 L.tprod3oilf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto pro3t;
ivreg2 ln_oilf14 L.ln_oilf14 L.tprod3oilf14  
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto pro4t;
ivreg2 ln_oilf14 L.ln_oilf14 L.tprod3oilf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prm_gas L.prm_oil L.pr2_gas L.pr2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.prm_gas L.prm_oil L.pr2_gas L.pr2_oil) bw(auto) robust;
est sto pro5t;	
	
# delimit; 
set more off;
* Permian gas total initial;
ivreg2 ln_gas_captf14 L.ln_gas_captf14 L.tprod3gasf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3gasfNpr L.prod3oilf21 ), 
	first savefirst savefprefix(fs_prgt) endog(lphh_ce lpwti_ce) orthog(L.prod3gasfNpr L.prod3oilf21 ) bw(auto) robust;
est sto prgpt;
ivreg2 ln_gas_captf14 L.ln_gas_captf14 L.tprod3gasf14 lphh_ce lpwti_ce, bw(auto) robust;
est sto prglst; 
ivreg2 ln_gas_captf14 L.ln_gas_captf14 L.tprod3gasf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc ), 
	endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto prg1t;
ivreg2 ln_gas_captf14 L.ln_gas_captf14 L.tprod3gasf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prod3oilf21 L.prod3oilfNpr L.prod3gasf21 L.prod3gasfNpr ), 
	endog(lphh_ce lpwti_ce) orthog(L.prod3oilf21 L.prod3oilfNpr L.prod3gasf21 L.prod3gasfNpr ) bw(auto) robust;
est sto prg2t;
ivreg2 ln_gas_captf14 L.ln_gas_captf14 L.tprod3gasf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.ln2_gas_pftb21 L.ln2_oil_pftb21), 
	endog(lphh_ce lpwti_ce) orthog(L.ln2_gas_pftb21 L.ln2_oil_pftb21) bw(auto) robust;
est sto prg3t;
ivreg2 ln_gas_captf14 L.ln_gas_captf14 L.tprod3gasf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21), 
	endog(lphh_ce lpwti_ce) orthog(L.totgas21 L.totoil21 L.totgaswell21 L.totoilwell21) bw(auto) robust;
est sto prg4t;
ivreg2 ln_gas_captf14 L.ln_gas_captf14 L.tprod3gasf14 
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.prm_gas L.prm_oil L.pr2_gas L.pr2_oil), 
	endog(lphh_ce lpwti_ce) orthog(L.prm_gas L.prm_oil L.pr2_gas L.pr2_oil) bw(auto) robust;
est sto prg5t;



* Now estimate sum of our 5 basins and all basins in DI;

* All five of our basins aggregated;
# delimit;
ivreg2 oilfot L.oilfot L.tprod3oilfot
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc L.tprod3oilf21 L.tprod3gasf21 ),
	first savefirst savefprefix(fs_5to) endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto otoiltp;
ivreg2 gasfot L.gasfot L.tprod3gasfot
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsh bsc  L.tprod3gasf21),
	first savefirst savefprefix(fs_5tg) endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto otgastp;

* sum of all basins in DI; 
ivreg2 oilftot L.oilftot L.tprod3oilftot
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent blhdd blcdd bsc bsh ),
	first savefirst savefprefix(fs_allto) endog(lphh_ce lpwti_ce) orthog() bw(auto) robust;
est sto alloiltp;
ivreg2 gasftot L.gasftot L.tprod3gasftot
	(lphh_ce lpwti_ce = lhdd_dev L.sum12_hdd_dev lcdd_dev L.sum12_cdd_dev lbrent),
	first savefirst savefprefix(fs_alltg) endog(lphh_ce lpwti_ce) bw(auto) robust;
est sto allgastp;



*************** Table 6 - Target initial production, preferred results;
# delimit;
esttab anop chop etop fwop prop otoilp alloilp , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp , labels( "KP F-stat" "Hansen's J" )) 
							rename(L.ln_oil_targf1 L.y L.ln_gas_targf1 L.y 
							L.ln_oil_targf5 L.y L.ln_gas_targf5 L.y
							L.ln_oil_targf8 L.y L.ln_gas_targf8 L.y
							L.ln_oil_targf9 L.y L.ln_gas_targf9 L.y
							L.ln_oil_targf14 L.y L.ln_gas_targf14 L.y
							L.oil_targfot L.y L.gas_targfot L.y
							L.oil_targftot L.y L.gas_targftot L.y
							L.prod3oilf1 L.prod L.prod3gasf1 L.prod
							L.prod3oilf5 L.prod L.prod3gasf5 L.prod
							L.prod3oilf8 L.prod L.prod3gasf8 L.prod
							L.prod3oilf9 L.prod L.prod3gasf9 L.prod
							L.prod3oilf14 L.prod L.prod3gasf14 L.prod
							L.prod3oilfot L.prod L.prod3gasfot L.prod
							L.prod3oilftot L.prod L.prod3gasftot L.prod)
							title({\b Target initial production:} {\i Two stage least squares estimates})
					mtitles("Anadarko" "Chautauqua" "East Texas" "Fort Worth" "Permian" "All 5" "All Basins")
					mgroups("Crude oil", pattern(1 0 0 0 0 0 0));
esttab angp chgp etgp fwgp prgp otgasp allgasp ,
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp , labels( "KP F-stat" "Hansen's J" )) 
						rename(L.ln_oil_targf1 L.y L.ln_gas_targf1 L.y 
							L.ln_oil_targf5 L.y L.ln_gas_targf5 L.y
							L.ln_oil_targf8 L.y L.ln_gas_targf8 L.y
							L.ln_oil_targf9 L.y L.ln_gas_targf9 L.y
							L.ln_oil_targf14 L.y L.ln_gas_targf14 L.y
							L.oil_targfot L.y L.gas_targfot L.y
							L.oil_targftot L.y L.gas_targftot L.y
							L.prod3oilf1 L.prod L.prod3gasf1 L.prod
							L.prod3oilf5 L.prod L.prod3gasf5 L.prod
							L.prod3oilf8 L.prod L.prod3gasf8 L.prod
							L.prod3oilf9 L.prod L.prod3gasf9 L.prod
							L.prod3oilf14 L.prod L.prod3gasf14 L.prod
							L.prod3oilfot L.prod L.prod3gasfot L.prod
							L.prod3oilftot L.prod L.prod3gasftot L.prod)
					mtitles("Anadarko" "Chautauqua" "East Texas" "Fort Worth" "Permian" "All 5" "All Basins")
					mgroups("Natural gas", pattern(1 0 0 0 0 0 0));
					


*************** Table 7 - Total initial production, preferred results;
# delimit;
esttab anopt chopt etopt fwopt propt otoiltp alloiltp , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat(idp widstat jp cstatp, labels("Underid" "WeakID" "J" "C stat")) 
							rename(L.ln_oilf1 L.y L.ln_gas_captf1 L.y 
							L.ln_oilf5 L.y L.ln_gas_captf5 L.y
							L.ln_oilf8 L.y L.ln_gas_captf8 L.y
							L.ln_oilf9 L.y L.ln_gas_captf9 L.y
							L.ln_oilf14 L.y L.ln_gas_captf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.tprod3oilf1 L.prod L.tprod3gasf1 L.prod
							L.tprod3oilf5 L.prod L.tprod3gasf5 L.prod
							L.tprod3oilf8 L.prod L.tprod3gasf8 L.prod
							L.tprod3oilf9 L.prod L.tprod3gasf9 L.prod
							L.tprod3oilf14 L.prod L.tprod3gasf14 L.prod
							L.tprod3oilfot L.prod L.tprod3gasfot L.prod
							L.tprod3oilftot L.prod L.tprod3gasftot L.prod)
							title({\b Total initial production of crude oil:} {\i 2SLS Estimates})
					mtitles("Anadarko" "Chautauqua" "East Texas" "Fort Worth" "Permian" "All 5" "All Basins");
esttab angpt chgpt etgpt fwgpt prgpt otgastp allgastp ,
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat(idp widstat jp cstatp, labels("Underid" "WeakID" "J" "C stat")) 
						rename(L.ln_oilf1 L.y L.ln_gas_captf1 L.y 
							L.ln_oilf5 L.y L.ln_gas_captf5 L.y
							L.ln_oilf8 L.y L.ln_gas_captf8 L.y
							L.ln_oilf9 L.y L.ln_gas_captf9 L.y
							L.ln_oilf14 L.y L.ln_gas_captf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.tprod3oilf1 L.prod L.tprod3gasf1 L.prod
							L.tprod3oilf5 L.prod L.tprod3gasf5 L.prod
							L.tprod3oilf8 L.prod L.tprod3gasf8 L.prod
							L.tprod3oilf9 L.prod L.tprod3gasf9 L.prod
							L.tprod3oilf14 L.prod L.tprod3gasf14 L.prod
							L.tprod3oilfot L.prod L.tprod3gasfot L.prod
							L.tprod3oilftot L.prod L.tprod3gasftot L.prod)
							title({\b Total initial production of natural gas:} {\i 2SLS Estimates})
					mtitles("Anadarko" "Chautauqua" "East Texas" "Fort Worth" "Permian" "All 5" "All Basins");



*************** Appendix Tables; 

******** First stage results from preferred specifications, Appendix Tables A1 to A5;
# delimit; 
* Anadarko; 
esttab fs_anolpwti_ce fs_anolphh_ce fs_anglpwti_ce fs_anglphh_ce 
		fs_anotlpwti_ce fs_anotlphh_ce fs_angtlpwti_ce fs_angtlphh_ce,
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress se
						rename(L.ln_oil_targf1 L.y L.ln_gas_targf1 L.y 
							L.ln_oilf1 L.y L.ln_gas_captf1 L.y
							L.prod3oilf1 L.prod L.prod3gasf1 L.prod
							L.tprod3oilf1 L.prod L.tprod3gasf1 L.prod)
					title({\b First stage coefficients:} {\i Anadarko Basin})
					mtitles("WTI" "HH" "WTI" "HH" "WTI" "HH" "WTI" "HH")
					mgroups("Target Oil" "Target Gas" "Total Oil" "Total Gas", pattern(1 0 1 0 1 0 1 0));
# delimit; 					
* Chautauqua; 
esttab fs_cholpwti_ce fs_cholphh_ce fs_chglpwti_ce fs_chglphh_ce 
		fs_chotlpwti_ce fs_chotlphh_ce fs_chgtlpwti_ce fs_chgtlphh_ce,
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress se 
						rename(L.ln_oil_targf5 L.y L.ln_gas_targf5 L.y 
							L.ln_oilf5 L.y L.ln_gas_captf5 L.y
							L.prod3oilf5 L.prod L.prod3gasf5 L.prod
							L.tprod3oilf5 L.prod L.tprod3gasf5 L.prod)
					title({\b First stage coefficients:} {\i Chautauqua Platform Basin})
					mtitles("WTI" "HH" "WTI" "HH" "WTI" "HH" "WTI" "HH")
					mgroups("Target Oil" "Target Gas" "Total Oil" "Total Gas", pattern(1 0 1 0 1 0 1 0));
# delimit; 
* East Texas; 
esttab fs_etolpwti_ce fs_etolphh_ce fs_etglpwti_ce fs_etglphh_ce 
		fs_etotlpwti_ce fs_etotlphh_ce fs_etgtlpwti_ce fs_etgtlphh_ce,
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress se 
						rename(L.ln_oil_targf8 L.y L.ln_gas_targf8 L.y 
							L.ln_oilf8 L.y L.ln_gas_captf8 L.y
							L.prod3oilf8 L.prod L.prod3gasf8 L.prod
							L.tprod3oilf8 L.prod L.tprod3gasf8 L.prod)
					title({\b First stage coefficients:} {\i East Texas Basin})
					mtitles("WTI" "HH" "WTI" "HH" "WTI" "HH" "WTI" "HH")
					mgroups("Target Oil" "Target Gas" "Total Oil" "Total Gas", pattern(1 0 1 0 1 0 1 0));
# delimit; 					
* Fort Worth; 
esttab fs_fwolpwti_ce fs_fwolphh_ce fs_fwglpwti_ce fs_fwglphh_ce 
		fs_fwotlpwti_ce fs_fwotlphh_ce fs_fwgtlpwti_ce fs_fwgtlphh_ce,
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress se 
						rename(L.ln_oil_targf9 L.y L.ln_gas_targf9 L.y 
							L.ln_oilf9 L.y L.ln_gas_captf9 L.y
							L.prod3oilf9 L.prod L.prod3gasf9 L.prod
							L.tprod3oilf9 L.prod L.tprod3gasf9 L.prod)
					title({\b First stage coefficients:} {\i Fort Worth Basin})
					mtitles("WTI" "HH" "WTI" "HH" "WTI" "HH" "WTI" "HH")
					mgroups("Target Oil" "Target Gas" "Total Oil" "Total Gas", pattern(1 0 1 0 1 0 1 0));
# delimit; 					
* Permian; 
esttab fs_prolpwti_ce fs_prolphh_ce fs_prglpwti_ce fs_prglphh_ce 
		fs_protlpwti_ce fs_protlphh_ce fs_prgtlpwti_ce fs_prgtlphh_ce,
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress se
						rename(L.ln_oil_targf14 L.y L.ln_gas_targf14 L.y 
							L.ln_oilf14 L.y L.ln_gas_captf14 L.y
							L.prod3oilf14 L.prod L.prod3gasf14 L.prod
							L.tprod3oilf14 L.prod L.tprod3gasf14 L.prod)
					title({\b First stage coefficients:} {\i Permian Basin})
					mtitles("WTI" "HH" "WTI" "HH" "WTI" "HH" "WTI" "HH")
					mgroups("Target Oil" "Target Gas" "Total Oil" "Total Gas", pattern(1 0 1 0 1 0 1 0));



******** Robustness to alternative IVs, Appendix Tables A6 to A15;
# delimit;
* Appendix Table A6 - Anadarko target initial, robustness to IVs;
esttab anop anols ano1 ano2 ano3 ano4 ano5 , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat")) 
							rename(L.ln_oil_targf1 L.y L.ln_gas_targf1 L.y 
							L.ln_oil_targf5 L.y L.ln_gas_targf5 L.y
							L.ln_oil_targf8 L.y L.ln_gas_targf8 L.y
							L.ln_oil_targf9 L.y L.ln_gas_targf9 L.y
							L.ln_oil_targf14 L.y L.ln_gas_targf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.prod3oilf1 L.prod L.prod3gasf1 L.prod
							L.prod3oilf5 L.prod L.prod3gasf5 L.prod
							L.prod3oilf8 L.prod L.prod3gasf8 L.prod
							L.prod3oilf9 L.prod L.prod3gasf9 L.prod
							L.prod3oilf14 L.prod L.prod3gasf14 L.prod
							L.prod3oilfot L.prod L.prod3gasfot L.prod
							L.prod3oilftot L.prod L.prod3gasftot L.prod)
							title({Anadarko target initial production of crude oil, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");
esttab angp angls ang1 ang2 ang3 ang4 ang5 , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat"))
							rename(L.ln_oil_targf1 L.y L.ln_gas_targf1 L.y 
							L.ln_oil_targf5 L.y L.ln_gas_targf5 L.y
							L.ln_oil_targf8 L.y L.ln_gas_targf8 L.y
							L.ln_oil_targf9 L.y L.ln_gas_targf9 L.y
							L.ln_oil_targf14 L.y L.ln_gas_targf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.prod3oilf1 L.prod L.prod3gasf1 L.prod
							L.prod3oilf5 L.prod L.prod3gasf5 L.prod
							L.prod3oilf8 L.prod L.prod3gasf8 L.prod
							L.prod3oilf9 L.prod L.prod3gasf9 L.prod
							L.prod3oilf14 L.prod L.prod3gasf14 L.prod
							L.prod3oilfot L.prod L.prod3gasfot L.prod
							L.prod3oilftot L.prod L.prod3gasftot L.prod)
							title({Anadarko target initial production of natural gas, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");

					

# delimit;
* Appendix Table A7 - Chautauqua Platform target initial, robustness to IVs;
esttab chop chols cho1 cho2 cho3 cho4 cho5, 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat")) 
							rename(L.ln_oil_targf1 L.y L.ln_gas_targf1 L.y 
							L.ln_oil_targf5 L.y L.ln_gas_targf5 L.y
							L.ln_oil_targf8 L.y L.ln_gas_targf8 L.y
							L.ln_oil_targf9 L.y L.ln_gas_targf9 L.y
							L.ln_oil_targf14 L.y L.ln_gas_targf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.prod3oilf1 L.prod L.prod3gasf1 L.prod
							L.prod3oilf5 L.prod L.prod3gasf5 L.prod
							L.prod3oilf8 L.prod L.prod3gasf8 L.prod
							L.prod3oilf9 L.prod L.prod3gasf9 L.prod
							L.prod3oilf14 L.prod L.prod3gasf14 L.prod
							L.prod3oilfot L.prod L.prod3gasfot L.prod
							L.prod3oilftot L.prod L.prod3gasftot L.prod)
							title({Chautauqua Platform target initial production of crude oil, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");
esttab chgp chgls chg1 chg2 chg3 chg4 chg5 , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat"))
							rename(L.ln_oil_targf1 L.y L.ln_gas_targf1 L.y 
							L.ln_oil_targf5 L.y L.ln_gas_targf5 L.y
							L.ln_oil_targf8 L.y L.ln_gas_targf8 L.y
							L.ln_oil_targf9 L.y L.ln_gas_targf9 L.y
							L.ln_oil_targf14 L.y L.ln_gas_targf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.prod3oilf1 L.prod L.prod3gasf1 L.prod
							L.prod3oilf5 L.prod L.prod3gasf5 L.prod
							L.prod3oilf8 L.prod L.prod3gasf8 L.prod
							L.prod3oilf9 L.prod L.prod3gasf9 L.prod
							L.prod3oilf14 L.prod L.prod3gasf14 L.prod
							L.prod3oilfot L.prod L.prod3gasfot L.prod
							L.prod3oilftot L.prod L.prod3gasftot L.prod)
							title({Chautauqua Platform target initial production of natural gas, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");


					
# delimit;
* Appendix Table A8 - East Texas target initial, robustness to IVs;
esttab etop etols eto1 eto2 eto3 eto4 eto5 , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat")) 
							rename(L.ln_oil_targf1 L.y L.ln_gas_targf1 L.y 
							L.ln_oil_targf5 L.y L.ln_gas_targf5 L.y
							L.ln_oil_targf8 L.y L.ln_gas_targf8 L.y
							L.ln_oil_targf9 L.y L.ln_gas_targf9 L.y
							L.ln_oil_targf14 L.y L.ln_gas_targf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.prod3oilf1 L.prod L.prod3gasf1 L.prod
							L.prod3oilf5 L.prod L.prod3gasf5 L.prod
							L.prod3oilf8 L.prod L.prod3gasf8 L.prod
							L.prod3oilf9 L.prod L.prod3gasf9 L.prod
							L.prod3oilf14 L.prod L.prod3gasf14 L.prod
							L.prod3oilfot L.prod L.prod3gasfot L.prod
							L.prod3oilftot L.prod L.prod3gasftot L.prod)
							title({East Texas target initial production of crude oil, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");
esttab etgp etgls etg1 etg2 etg3 etg4 etg5 , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat")) 
							rename(L.ln_oil_targf1 L.y L.ln_gas_targf1 L.y 
							L.ln_oil_targf5 L.y L.ln_gas_targf5 L.y
							L.ln_oil_targf8 L.y L.ln_gas_targf8 L.y
							L.ln_oil_targf9 L.y L.ln_gas_targf9 L.y
							L.ln_oil_targf14 L.y L.ln_gas_targf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.prod3oilf1 L.prod L.prod3gasf1 L.prod
							L.prod3oilf5 L.prod L.prod3gasf5 L.prod
							L.prod3oilf8 L.prod L.prod3gasf8 L.prod
							L.prod3oilf9 L.prod L.prod3gasf9 L.prod
							L.prod3oilf14 L.prod L.prod3gasf14 L.prod
							L.prod3oilfot L.prod L.prod3gasfot L.prod
							L.prod3oilftot L.prod L.prod3gasftot L.prod)
							title({East Texas target initial production of natural gas, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");

					

# delimit;
* Appendix Table A9 - Fort Worth target initial, robustness to IVs;
esttab fwop fwols fwo1 fwo2 fwo3 fwo4 fwo5 , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat"))
							rename(L.ln_oil_targf1 L.y L.ln_gas_targf1 L.y 
							L.ln_oil_targf5 L.y L.ln_gas_targf5 L.y
							L.ln_oil_targf8 L.y L.ln_gas_targf8 L.y
							L.ln_oil_targf9 L.y L.ln_gas_targf9 L.y
							L.ln_oil_targf14 L.y L.ln_gas_targf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.prod3oilf1 L.prod L.prod3gasf1 L.prod
							L.prod3oilf5 L.prod L.prod3gasf5 L.prod
							L.prod3oilf8 L.prod L.prod3gasf8 L.prod
							L.prod3oilf9 L.prod L.prod3gasf9 L.prod
							L.prod3oilf14 L.prod L.prod3gasf14 L.prod
							L.prod3oilfot L.prod L.prod3gasfot L.prod
							L.prod3oilftot L.prod L.prod3gasftot L.prod)
							title({Fort Worth target initial production of crude oil, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");
esttab fwgp fwgls fwg1 fwg2 fwg3 fwg4 fwg5 , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat")) 
							rename(L.ln_oil_targf1 L.y L.ln_gas_targf1 L.y 
							L.ln_oil_targf5 L.y L.ln_gas_targf5 L.y
							L.ln_oil_targf8 L.y L.ln_gas_targf8 L.y
							L.ln_oil_targf9 L.y L.ln_gas_targf9 L.y
							L.ln_oil_targf14 L.y L.ln_gas_targf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.prod3oilf1 L.prod L.prod3gasf1 L.prod
							L.prod3oilf5 L.prod L.prod3gasf5 L.prod
							L.prod3oilf8 L.prod L.prod3gasf8 L.prod
							L.prod3oilf9 L.prod L.prod3gasf9 L.prod
							L.prod3oilf14 L.prod L.prod3gasf14 L.prod
							L.prod3oilfot L.prod L.prod3gasfot L.prod
							L.prod3oilftot L.prod L.prod3gasftot L.prod)
							title({Fort Worth target initial production of natural gas, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");
					


# delimit;
* Appendix Table A10 - Permian target initial, robustness to IVs;
esttab prop prols pro1 pro2 pro3 pro4 pro5 , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat")) 
							rename(L.ln_oil_targf1 L.y L.ln_gas_targf1 L.y 
							L.ln_oil_targf5 L.y L.ln_gas_targf5 L.y
							L.ln_oil_targf8 L.y L.ln_gas_targf8 L.y
							L.ln_oil_targf9 L.y L.ln_gas_targf9 L.y
							L.ln_oil_targf14 L.y L.ln_gas_targf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.prod3oilf1 L.prod L.prod3gasf1 L.prod
							L.prod3oilf5 L.prod L.prod3gasf5 L.prod
							L.prod3oilf8 L.prod L.prod3gasf8 L.prod
							L.prod3oilf9 L.prod L.prod3gasf9 L.prod
							L.prod3oilf14 L.prod L.prod3gasf14 L.prod
							L.prod3oilfot L.prod L.prod3gasfot L.prod
							L.prod3oilftot L.prod L.prod3gasftot L.prod)
							title({Permian target initial production of crude oil, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");
esttab prgp prgls prg1 prg2 prg3 prg4 prg5 , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat")) 
							rename(L.ln_oil_targf1 L.y L.ln_gas_targf1 L.y 
							L.ln_oil_targf5 L.y L.ln_gas_targf5 L.y
							L.ln_oil_targf8 L.y L.ln_gas_targf8 L.y
							L.ln_oil_targf9 L.y L.ln_gas_targf9 L.y
							L.ln_oil_targf14 L.y L.ln_gas_targf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.prod3oilf1 L.prod L.prod3gasf1 L.prod
							L.prod3oilf5 L.prod L.prod3gasf5 L.prod
							L.prod3oilf8 L.prod L.prod3gasf8 L.prod
							L.prod3oilf9 L.prod L.prod3gasf9 L.prod
							L.prod3oilf14 L.prod L.prod3gasf14 L.prod
							L.prod3oilfot L.prod L.prod3gasfot L.prod
							L.prod3oilftot L.prod L.prod3gasftot L.prod)
							title({Permian target initial production of natural gas, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");
			
			

# delimit;
* Appendix Table A11 - Anadarko total initial, robustness to IVs;
esttab anopt anolst ano1t ano2t ano3t ano4t ano5t , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat")) 
							rename(L.ln_oilf1 L.y L.ln_gas_captf1 L.y 
							L.ln_oilf5 L.y L.ln_gas_captf5 L.y
							L.ln_oilf8 L.y L.ln_gas_captf8 L.y
							L.ln_oilf9 L.y L.ln_gas_captf9 L.y
							L.ln_oilf14 L.y L.ln_gas_captf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.tprod3oilf1 L.prod L.tprod3gasf1 L.prod
							L.tprod3oilf5 L.prod L.tprod3gasf5 L.prod
							L.tprod3oilf8 L.prod L.tprod3gasf8 L.prod
							L.tprod3oilf9 L.prod L.tprod3gasf9 L.prod
							L.tprod3oilf14 L.prod L.tprod3gasf14 L.prod
							L.tprod3oilfot L.prod L.tprod3gasfot L.prod
							L.tprod3oilftot L.prod L.tprod3gasftot L.prod)
							title({Anadarko total initial production of crude oil, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");
esttab angpt anglst ang1t ang2t ang3t ang4t ang5t , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat"))
							rename(L.ln_oilf1 L.y L.ln_gas_captf1 L.y 
							L.ln_oilf5 L.y L.ln_gas_captf5 L.y
							L.ln_oilf8 L.y L.ln_gas_captf8 L.y
							L.ln_oilf9 L.y L.ln_gas_captf9 L.y
							L.ln_oilf14 L.y L.ln_gas_captf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.tprod3oilf1 L.prod L.tprod3gasf1 L.prod
							L.tprod3oilf5 L.prod L.tprod3gasf5 L.prod
							L.tprod3oilf8 L.prod L.tprod3gasf8 L.prod
							L.tprod3oilf9 L.prod L.tprod3gasf9 L.prod
							L.tprod3oilf14 L.prod L.tprod3gasf14 L.prod
							L.tprod3oilfot L.prod L.tprod3gasfot L.prod
							L.tprod3oilftot L.prod L.tprod3gasftot L.prod)
							title({Anadarko total initial production of natural gas, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");

					
					
# delimit;
* Appendix Table A12 - Chautauqua Platform total initial, robustness to IVs;
esttab chopt cholst cho1t cho2t cho3t cho4t cho5t , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat")) 
							rename(L.ln_oilf1 L.y L.ln_gas_captf1 L.y 
							L.ln_oilf5 L.y L.ln_gas_captf5 L.y
							L.ln_oilf8 L.y L.ln_gas_captf8 L.y
							L.ln_oilf9 L.y L.ln_gas_captf9 L.y
							L.ln_oilf14 L.y L.ln_gas_captf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.tprod3oilf1 L.prod L.tprod3gasf1 L.prod
							L.tprod3oilf5 L.prod L.tprod3gasf5 L.prod
							L.tprod3oilf8 L.prod L.tprod3gasf8 L.prod
							L.tprod3oilf9 L.prod L.tprod3gasf9 L.prod
							L.tprod3oilf14 L.prod L.tprod3gasf14 L.prod
							L.tprod3oilfot L.prod L.tprod3gasfot L.prod
							L.tprod3oilftot L.prod L.tprod3gasftot L.prod)
							title({Chautauqua Platform total initial production of crude oil, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");
esttab chgpt chglst chg1t chg2t chg3t chg4t chg5t , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat"))
							rename(L.ln_oilf1 L.y L.ln_gas_captf1 L.y 
							L.ln_oilf5 L.y L.ln_gas_captf5 L.y
							L.ln_oilf8 L.y L.ln_gas_captf8 L.y
							L.ln_oilf9 L.y L.ln_gas_captf9 L.y
							L.ln_oilf14 L.y L.ln_gas_captf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.tprod3oilf1 L.prod L.tprod3gasf1 L.prod
							L.tprod3oilf5 L.prod L.tprod3gasf5 L.prod
							L.tprod3oilf8 L.prod L.tprod3gasf8 L.prod
							L.tprod3oilf9 L.prod L.tprod3gasf9 L.prod
							L.tprod3oilf14 L.prod L.tprod3gasf14 L.prod
							L.tprod3oilfot L.prod L.tprod3gasfot L.prod
							L.tprod3oilftot L.prod L.tprod3gasftot L.prod)
							title({Chautauqua Platform total initial production of natural gas, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");

					
					
# delimit;
* Appendix Table A13 - East Texas total initial, robustness to IVs;
esttab etopt etolst eto1t eto2t eto3t eto4t eto5t, 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat")) 
							rename(L.ln_oilf1 L.y L.ln_gas_captf1 L.y 
							L.ln_oilf5 L.y L.ln_gas_captf5 L.y
							L.ln_oilf8 L.y L.ln_gas_captf8 L.y
							L.ln_oilf9 L.y L.ln_gas_captf9 L.y
							L.ln_oilf14 L.y L.ln_gas_captf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.tprod3oilf1 L.prod L.tprod3gasf1 L.prod
							L.tprod3oilf5 L.prod L.tprod3gasf5 L.prod
							L.tprod3oilf8 L.prod L.tprod3gasf8 L.prod
							L.tprod3oilf9 L.prod L.tprod3gasf9 L.prod
							L.tprod3oilf14 L.prod L.tprod3gasf14 L.prod
							L.tprod3oilfot L.prod L.tprod3gasfot L.prod
							L.tprod3oilftot L.prod L.tprod3gasftot L.prod)
							title({East Texas total initial production of crude oil, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");
esttab etgpt etglst etg1t etg2t etg3t etg4t etg5t , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat")) 
							rename(L.ln_oilf1 L.y L.ln_gas_captf1 L.y 
							L.ln_oilf5 L.y L.ln_gas_captf5 L.y
							L.ln_oilf8 L.y L.ln_gas_captf8 L.y
							L.ln_oilf9 L.y L.ln_gas_captf9 L.y
							L.ln_oilf14 L.y L.ln_gas_captf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.tprod3oilf1 L.prod L.tprod3gasf1 L.prod
							L.tprod3oilf5 L.prod L.tprod3gasf5 L.prod
							L.tprod3oilf8 L.prod L.tprod3gasf8 L.prod
							L.tprod3oilf9 L.prod L.tprod3gasf9 L.prod
							L.tprod3oilf14 L.prod L.tprod3gasf14 L.prod
							L.tprod3oilfot L.prod L.tprod3gasfot L.prod
							L.tprod3oilftot L.prod L.tprod3gasftot L.prod)
							title({East Texas total initial production of natural gas, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");
					
					
					
# delimit;
* Appendix Table A14 - Fort Worth total initial, robustness to IVs;
esttab fwopt fwolst fwo1t fwo2t fwo3t fwo4t fwo5t , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat"))
							rename(L.ln_oilf1 L.y L.ln_gas_captf1 L.y 
							L.ln_oilf5 L.y L.ln_gas_captf5 L.y
							L.ln_oilf8 L.y L.ln_gas_captf8 L.y
							L.ln_oilf9 L.y L.ln_gas_captf9 L.y
							L.ln_oilf14 L.y L.ln_gas_captf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.tprod3oilf1 L.prod L.tprod3gasf1 L.prod
							L.tprod3oilf5 L.prod L.tprod3gasf5 L.prod
							L.tprod3oilf8 L.prod L.tprod3gasf8 L.prod
							L.tprod3oilf9 L.prod L.tprod3gasf9 L.prod
							L.tprod3oilf14 L.prod L.tprod3gasf14 L.prod
							L.tprod3oilfot L.prod L.tprod3gasfot L.prod
							L.tprod3oilftot L.prod L.tprod3gasftot L.prod)
							title({Fort Worth total initial production of crude oil, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");
esttab fwgpt fwglst fwg1t fwg2t fwg3t fwg4t fwg5t , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat")) 
							rename(L.ln_oilf1 L.y L.ln_gas_captf1 L.y 
							L.ln_oilf5 L.y L.ln_gas_captf5 L.y
							L.ln_oilf8 L.y L.ln_gas_captf8 L.y
							L.ln_oilf9 L.y L.ln_gas_captf9 L.y
							L.ln_oilf14 L.y L.ln_gas_captf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.tprod3oilf1 L.prod L.tprod3gasf1 L.prod
							L.tprod3oilf5 L.prod L.tprod3gasf5 L.prod
							L.tprod3oilf8 L.prod L.tprod3gasf8 L.prod
							L.tprod3oilf9 L.prod L.tprod3gasf9 L.prod
							L.tprod3oilf14 L.prod L.tprod3gasf14 L.prod
							L.tprod3oilfot L.prod L.tprod3gasfot L.prod
							L.tprod3oilftot L.prod L.tprod3gasftot L.prod)
							title({Fort Worth total initial production of natural gas, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");
		
		
		
# delimit;
* Appendix Table A15 - Permian total initial, robustness to IVs;
esttab propt prolst pro1t pro2t pro3t pro4t pro5t , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat")) 
							rename(L.ln_oilf1 L.y L.ln_gas_captf1 L.y 
							L.ln_oilf5 L.y L.ln_gas_captf5 L.y
							L.ln_oilf8 L.y L.ln_gas_captf8 L.y
							L.ln_oilf9 L.y L.ln_gas_captf9 L.y
							L.ln_oilf14 L.y L.ln_gas_captf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.tprod3oilf1 L.prod L.tprod3gasf1 L.prod
							L.tprod3oilf5 L.prod L.tprod3gasf5 L.prod
							L.tprod3oilf8 L.prod L.tprod3gasf8 L.prod
							L.tprod3oilf9 L.prod L.tprod3gasf9 L.prod
							L.tprod3oilf14 L.prod L.tprod3gasf14 L.prod
							L.tprod3oilfot L.prod L.tprod3gasfot L.prod
							L.tprod3oilftot L.prod L.tprod3gasftot L.prod)
							title({Permian total initial production of crude oil, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");		
esttab prgpt prglst prg1t prg2t prg3t prg4t prg5t , 
					nonumber star(* 0.1 ** 0.05 *** 0.01) nogaps 
					b(a2) se(a2) compress r2 se 
					stat( widstat jp cstatp, labels("KP F-stat" "Hansen's J" "C stat")) 
							rename(L.ln_oilf1 L.y L.ln_gas_captf1 L.y 
							L.ln_oilf5 L.y L.ln_gas_captf5 L.y
							L.ln_oilf8 L.y L.ln_gas_captf8 L.y
							L.ln_oilf9 L.y L.ln_gas_captf9 L.y
							L.ln_oilf14 L.y L.ln_gas_captf14 L.y
							L.oilfot L.y L.gasfot L.y
							L.oilftot L.y L.gasftot L.y
							L.tprod3oilf1 L.prod L.tprod3gasf1 L.prod
							L.tprod3oilf5 L.prod L.tprod3gasf5 L.prod
							L.tprod3oilf8 L.prod L.tprod3gasf8 L.prod
							L.tprod3oilf9 L.prod L.tprod3gasf9 L.prod
							L.tprod3oilf14 L.prod L.tprod3gasf14 L.prod
							L.tprod3oilfot L.prod L.tprod3gasfot L.prod
							L.tprod3oilftot L.prod L.tprod3gasftot L.prod)
							title({Permian total initial production of natural gas, Comparing 2SLS Estimates})
					mtitles("Preferred" "OLS" "Weather/Brent" "Q/Well" "Q/Depth" "Cumulative" "Non-Neighbor Q");



